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Abstract 



We investigate the usefulness of the "string-inspired technique" for gauge the- 
ory calculations in a constant external field background. Our approach is based 
on Strassler's worldline path integral approach to the Bern-Kosower formalism, 
and on the construction of worldline (super-) Green's functions incorporating 
external fields as well as internal propagators. The worldline path integral rep- 
resentation of the gluon loop is reexamined in detail. We calculate the two-loop 
effective actions induced for a constant external field by a scalar and spinor 
loop, and the corresponding one-loop effective action in the gluon loop case. 
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1. Introduction 



It is by now well-established that techniques from string perturbation theory can be 
used to improve on calculational efficiency in ordinary quantum field theory. The rele- 
vance of string theory for this purpose is based on the fact that many, and perhaps all, 
amplitudes in quantum field theory can be represented as the infinite string tension limits 
of appropriately chosen (super-) string amplitudes. This is, of course, an intrinsic and 
well-known property of string-theory. 

It is a more recent discovery, however, that such representations can lead to an in- 
teresting alternative to standard Feynman diagram calculations. Following earlier work 
on the P - function for Yang-Mills-theory ^, ^, Bern and Kosower Q considered 
the infinite string tension limit of gauge boson scattering amplitudes, formulated in an 
appropriate heterotic string model. By a detailed analysis of this limit, they succeeded in 
deriving a novel type of "Feynman Rules" for the construction of ordinary one-loop gluon 
(photon) scattering amplitudes. The resulting integral representations are equivalent to 
the ones originating from Feynman diagrams |^] , but lead to a significant reduction of the 
number of terms to be computed in gauge theory calculations. This property was then 
successfully exploited to obtain both five gluon and four graviton amplitudes 

Strassler |^ later showed that, for many cases of interest, the same integral represen- 
tations can be derived from a first-quantized reformulation of ordinary field theory. In this 
approach, one starts with writing the one-loop effective action as a particle (super) path- 
integral, of a type which has been known for many years [|lO|, |ll|, 12, 13, |l^, |l5|, |l^, p^, p^. 
Those path-integrals are then considered as the field-theory analogues of the Polyakov path 
integral, and evaluated in a way analogous to string perturbation theory (some suggestions 
along similar lines had already been made in p^ ). 

The resulting formalism offers an alternative to standard field theory techniques which 
circumvents much of the apparatus of quantum field theory. It works equally well for 
effective action and scattering amplitude calculations, on- and off-shell. It has been 
applied to a number of calculations in gauge theory |2^, 21, 22, |2^, 24] and generalized 
to cases where, besides gauge and scalar self-couplings, Yukawa [25, 26, 27 1 and axial 
couplings pr 



|, |2^, 27 1 are present. 

Due to its simplicity, it appears also well-suited to the construction of multiloop gener- 
alizations of the Bern-Kosower formalism. Steps in this direction have already been taken 
by various authors, and along different lines |^, |3^, ^ . In particular, the original Bern- 
Kosower program becomes very hard to carry through at the two-loop level, due to the 
complicated structure of genus two string amplitudes. Nevertheless, recently substantial 
progress has beeen achieved in this line of work 



|32 



A multiloop generalization following the spirit of Strassler's approach has been pro- 
posed by two of the present authors, first for scalar field theories in |34|. This general- 
ization uses the concept of worldline Green's functions on graphs [^, 35, 33], and leads 



to integral representations combining whole classes of graphs. This work was extended to 
QED in ||3^, and its practical viability demonstrated by a recalculation of the two- loop 
(scalar and spinor) QED (3 - functions. For the case of multiloop amplitudes in scalar 
QED, a more comprehensive treatment along the same lines was given in |37|. This 
includes amplitudes involving external scalars. 

An important role in quantum electrodynamics is played by calculations involving con- 
stant external fields. This subject originates with Euler and Heisenberg's classic one-loop 
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calculation of the static limit of photon scattering in spinor quantum electrodynamics 
[38|. Schwinger's introduction of the proper-time method in 1951 [39| allowed him to 
reproduce this result, and the analogous one for scalar quantum electrodynamics, with 
considerably less effort. For calculations of this type, it has turned out generally advanta- 
geous to take account of the external field already at the level of the Feynman rules, i.e. 
to absorb it into the free propagators. Appropriate formalisms have been developed both 



for QED |4C, 41, m and QCD |43]. 



The purpose of the present paper is threefold. First, we would like to cast the worldline 
multiloop formalism proposed in |34, 36| into a form suitable for constant external field 
calculations in quantum electrodynamics. This will be achieved in a way analogous to field 
theory, namely by modifying the worldline Green's functions so as to take the external 
field into account. At the one- loop level, similar proposals have already been made by 
several authors ||2|, |4|, Q (see also Q). 

Second, we will apply this formalism to a two-loop effective action calculation in QED. 
In [^] the worldline formalism was used for a recalculation of the Euler-Heisenberg 
Lagrangian. In the present work, we will extend this calculation to the two-loop level, 
and calculate the correction to this effective Lagrangian due to one internal photon, both 
for scalar and spinor QED. We discuss in how far this calculation improves on previous 



calculations of the same quantities M, 47, 45, 49 



Third, we reconsider the super path integral representation of the spin 1 particle, 
and use this path integral for calculating the effective action induced for an external 
constant pseudo-abelian field by a gluon loop in Yang-Mills- theory. While external gluons 
pose no particular problems in the worldline formalism, apart from forcing path-ordering, 
internal gluons are a much more delicate matter. While it is not difficult to construct 
free path integrals describing particles of arbitrary spin, those constructions usually run 
into consistency problems if one tries to couple a path integral with spin higher than ^ 
to a spin-1 background 52|. A somewhat non-standard path integral, mimicking the 
superstring, had been proposed by Strassler in his rederivation of the Bern-Kosower 
rules for the gluon loop case. To the best of our knowledge, it has not yet been proven that 
this path integral correctly reproduces all the pinch terms implicit in the Bern-Kosower 
master formula. 

We will first give a rigorous derivation of a spin-1 path integral which, while not 
identical with the one used by Strassler, is easily seen to be equivalent. We then use it 
for calculating the effective action induced by a gluon loop for a constant pseudo-abelian 
background gluon field. The result will again be in agreement with the literature |53], 
and provides a nontrivial check on the correctness of Strassler's proposal. 

The organization of this paper is as follows. In chapter 2 we review the worldline- 
formalism for one- loop photon scattering in scalar and spinor QED, and its generalization 
to gluon scattering. We then indicate the changes which are necessary to take a constant 
external field into account. Chapter 3 extends this analysis to the gluon loop. This so 
far includes worldline calculations of the one-loop effective actions induced for a (pseudo- 
abelian) constant external field by a scalar, spinor and gluon loop. For the QED case, we 
then extend the constant field formalism to the multiloop level in chapter 4. We apply it 
to calculations of the corresponding two-loop effective actions for scalar QED in chapter 
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5 and for spinor QED in chapter 6. Chapter 7 contains some remarks on the various 
ways of calculating the two-loop QED (3 - function in this formalism. Our results will be 
discussed in chapter 8. In appendix A we discuss the path integral representation of the 
electron propagator in an external field. In appendix B we derive the various worldline 
Green's functions used in our calculations, and discuss some of their properties. Appendix 
C contains a collection of some results concerning the determinants which we encounter 
in the evaluation of the spin-1 path integral. 



2. One Loop Photon Amplitudes in a Constant Background Field 

We begin with shortly reviewing how one-loop photon (gluon) scattering amplitudes 



are calculated in the worldline formalism 21, In his rederivation of the Bern- 
Kosower rules for gauge boson scattering off a spinor loop, Strassler sets out from the 
following well-known path integral representation for the corresponding one-loop effective 
action (see e.g. 54]): 



T[A] 



'dT 

/o 

xtvVexp 



Vx / Vip 



1 



(2.1) 



In this formula, T is the usual Schwinger proper-time parameter. The periodic functions 
x^{t) describe the embedding of the circle with circumference T into D - dimensional 
Euclidean spacetime, while the ^/^^(r)'s are antiperiodic Grassmann functions. The peri- 
odicity properties are expressed by the subscripts -P, A on the path integral. The colour 
trace tr and the path ordering V apply, of course, only to the nonabelian case. We have 
chosen a constant euclidean worldline metric. 

The case of a scalar loop is obtained simply by discarding all Grassmann quantities 
and the global factor of -2, which takes care of the difference in statistics and degrees of 
freedom. 

Analogous path integral representations exist for the scalar and electron propagators 

imin. 



in a background field |1C, 11, [S^]- The integration is then over a space of paths with 
fixed boundary conditions. Those are obvious in the scalar case, while in the fermionic 
case there has been some debate on the correct boundary conditions to be imposed on the 
Grassmann path integral 



55|, 56, 19 1 . An adequate path integral representation for the 
gluon propagator seems to be missing in the literature, and will be derived in chapter 3. 
The more familiar cases of the scalar and spin-^ propagators are included in appendix A 
for completeness. 

While in the present paper we will make use only of the loop path integrals, we ex- 
pect those propagator path integral representations to play an important role in future 
extensions of this formalism. 

A number of different techniques have been applied to calculating this worldline path 
integral, and various generalizations thereof (see e.g. [57,|5S, 59,^0, 61, 32,63, 64, |65|, |6^). 
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In the "string-inspired" approach, the first step is to split the coordinate path integral 
into center of mass and relative coordinates, 

J Vx = j dxQ J Vy 
/ dry^'ir) = . 

^0 

(2.2) 

The path integrals over y and ■0 are then evaluated by Wick contractions, as in a one- 
dimensional field theory on the circle. The Green's functions to be used are 



{y^{ri)y''{T2)) 

{ririW{r2)) 



-9' 



T 



^ g'"'GF{Tl,T2) = ^ 5^''sign(Ti - T2) 



(2.3) 



We will often abbreviate Gb{ti,T2) =: Gbi2 etc. They solve the differential equations 



1 02 



tGb{ti-,T2) 



^^GHn,r2) 



<5(ri - T2) - - 
5{ti - T2) 



(2.4) 



with periodic (antiperiodic) boundary conditions for Gb{Gf)- Some freedom exists in the 



definition of Gb, which has been discussed elsewhere ||22|, 67]; in particular, a constant 
added to would drop out after momentum conservation is applied. 
With our conventions, the free path integrals are normalized as 



A 



Vy exp 

p 

ViIj exp 



T I 
dT-f 

4 
T I . 

dr—iptp 
2 



(2.5) 



The result of this evaluation is the one-loop effective Lagrangian C{xo). Combined with a 
covariant Taylor expansion of the external field at xq , this yields a new and highly efficient 
algorithm for calculating higher derivative expansions in gauge theory |2^. 
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One-loop scattering amplitudes are obtained by specializing the background to a finite 
sum of plane waves of definite polarization. Equivalently, one may define integrated vertex 
operators 

Va = r dr [x^£^ - 2iiP^^^k^E^]e^Y>[ikx{T)] (2.6) 
Jo 

for external photons (gluons) of definite momentum and polarization, and calculate mul- 
tiple insertions of those vertex operators into the free path integral (of course, only the 
first term has to be taken for the scalar loop). denotes a generator of the gauge group 
in the representation of the loop particle. In the nonabelian case, for the spinor loop an 
additional two-gluon vertex operator is required Q. 

The path integrals are performed using the well-known formulas for Wick-contractions 
involving exponentials, e.g. 

exp[iA;i • x(ri)]exp[iA;2 • x{t2)]^ = exp GB{Ti,T2)ki ■ k2 (2.7) 

etc. (the factors x^e^ may be formally exponentiated for convenience). 

Explicit execution of the ip - path integral would be algebraically equivalent to the 
calculation of the corresponding Dirac traces in field theory. It can be circumvented by the 
following remarkable feature of the Bern-Kosower formalism, which may be understood 
as a consequence of worldsheet or worldline Q supersymmetry. After evaluation 
of the X - path integral, one is left with an integral over the parameters T,ti, . . . ,tn, 
where is the number of external legs. In the nonabelian case, the path-ordering leads 
to ordered r - integrals, / Hila^ dTi6{Ti — Tj+i). The integrand is an expression consisting 
of an exponential factor 

exp ^ Gb{ '^i 5 '^j ) ki kj 

i<j 

multiplied by a polynomial in the first and second derivatives of G^, 



GB{n,T2) = sign(ri - T2) - 2 

2 

G's(Ti,r2) = 2(5(ri - T2) - — 



(2.8) 



(here and in the following, a "dot" denotes differentiation with respect to the first variable) . 

All Gb's can be eliminated by partial integrations on the worldline, leading to an 
equivalent parameter integral dependent only on and Gb- According to the Bern- 
Kosower rules, all contributions from fermionic Wick contractions may then be taken into 
account simply by simultaneously replacing every closed cycle of G^'s appearing, say 
GBiiiiGBiih ■ ■ ■ GBinh, by its "supersymmetrization" , 

G Biii2G Bi2h ' ' ' Gsinh ~^ GBi^i2GBi2is ' ' ' GBinh ~ G Fiii2G Fi2i3 ' ' ' Gpinh- (2-9) 
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Note that an expression is considered a cycle already if it can be put into cycle form using 
the antisymmetry of Gb (e.g. GBabGBab = —GBabGBba)- Unfortunately the practical 
value of this procedure rapidly diminishes with increasing number of external legs (> 5), 
as the number of terms making up the integrand starts to significantly increase in the 
partial integration procedure |]69|| . 

The worldline supersymmetry (see eq. (A. 9)) makes it possible to combine the x - and 



Tp - path integrals into the following super path integral p3, 19, 55, 70 



T[A] 



■dT 



e "^^^ iVXiiVe 



■^XD'-^X-ieDX^,A^,(X) 



(2.10) 



D 



de dr 

1 



The photon (gluon) vertex operator is then rewritten as 



/ dTdeenDXr,e^^[ikX], 
Jo 



and we are left with a single Wick-contraction rule 



(2.11) 



(2.12) 



{Y''{n,ei)Y''{T2,02)) = -g^''Gin,ei;T2,e2) 

G{ri,ei;T2,e2) = GBiTl,T2) + e,e2GF{Ti,T2). 

(2.13) 

The fermion loop case can thus be made to look formally identical to the scalar loop case, 
and be regarded as its "supersymmetrization" . This analogy has its roots in the fact 
that the string-inspired technique corresponds to the use of a second-order formalism for 
fermions in field theory (see [^l| and ref. therein), instead of the usual first-order ones. In 
practical terms it means that the Grassmann Wick contractions are replaced by a number 
of Grassmann integrals, which have to be performed at a later stage of the calculation. 
Ultimately the superfield formalism leads to the same collection of parameter integrals 
to be performed, however we have found it useful for keeping intermediate expressions 
compact. In the nonabelian case it has the further advantage that the introduction of an 
additional two-gluon vertex operator can be avoided. Instead one introduces a suitable 
super symmetric generalization of the functions 0(rj — Tj+i) appearing in the ordered r - 
integrals [70|. 



Now let us assume that we have, in addition to the background field A^{x) we started 
with, a second one, A'^{x), with constant field strength tensor F^y. We will restrict 
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ourselves to the abelian case for the remainder of this chapter. Using Fock-Schwinger 
gauge centered at xq, we may take A^(x) to be of the form 

A^ix) = ^y,F,^ . (2.14) 
The constant field contribution to the worldline lagrangian may then be written as 

AC = ^ieyf^F^ytji, - ieipf.F^^ipy (2.15) 

in components, or as 

A£ = -h^eY^F^^DYu (2.16) 

in the superfield formalism. 

In any case, it is still quadratic in the worldline fields, and therefore need not be 
considered part of the interaction lagrangian; we can absorb it into the free worldline 
propagator (s). This means that we need to solve, instead of eqs. ( |2.4D for the worldline 
Green's functions, the modified equations 



1 



92 _ a 

2ieF— ]GB{n,T2) 



2 \dn^ 

^(^-^^-2ieF]gF{n,T2) 



<5(n - T2) 
Sin - T2) 



(2.17) 
(2.18) 



These equations will be solved in appendix B, with the result (deleting the "bar" again) 



aB(ri,r2) = ( . e-^^g^^^ + ieFGBi2 - ^ 

2{eF) \sm{eFT) T 



Q-ieFTGBV2 

Qf{ti,T2) = Gf12 



cos(eFr) 

(2.19) 



Equivalent expressions have been given for the pure magnetic field case in [p3| , and for the 
general case in p4|. These expressions should be understood as power series in the field 



strength matrix F (note that eqs.( |2.19 ) do not assume invertibility of F). Note also that 
the generalized Green's functions are still translation invariant in r, and thus functions of 
Ti — T2. They are, in general, nontrivial Lorentz matrices, so that the Wick contraction 



rules eq.(2.3) have to be rewritten as 



{y^iri)y''ir2)) = -ar(Ti,r2) 



1 



{r{n)r{r2)) = ^ar(n,r2). 



2" 

(2.20) 



7 



Again momentum conservation leads to the freedom to subtract from Qb its constant 
coincidence limit, 



SBiT,r) = ^^^^{eFcotieFT)--] (2.21) 



^-^f eFcot(eFr) - ^ 
To correctly obtain this and other coincidence limits, one has to apply the rules 



G'B(r,r) = 0, G|(r,r) = l. (2.22) 

More generally, coincidence limits should always be taken after derivatives. 
We also need the generalizations of Gb,Gb, which turn out to be 

Mr...) . 2(nl(<).-2„F,-|..) = A(_f£_,-~_i^ 

QB{r^.T2) ^ 2(ri|(I-2ieFa.-i) '|r2) = 25i2-2-— — e-^^™^^ 

sm(ei^ L ) 

(2.23) 

Let us also give the first few terms of the Taylor expansions in F for those four functions: 

Qbx2 = GBi2-^-'-GBi2GBi2TeF + {'^Gl,,-'^){eFf + 0{F^) 
^512 = Gbi2 + 2i{GBi2 - |)eF + ^GBi2GBi2T(eF)2 + 0{F'') 
Qbx2 = GBi2 + 2iGBi2eF-A{GBi2-'^){eFf + 0{F^) 





Qfi2 = GFi2-iGFi2GBi2TeF + 2GFi2GBi2T{eFY + 0{F'^ 



Again Qb and Qp may be assembled into a superpropagator. 



(2.24) 



a(ri,^i;T2,02) = GB{ri,T2)+9i92gF{Tl,T2). (2.25) 

At first sight this definition would seem not to accomodate the nonvanishing coincidence 
limit of ^iT' (which can not be subtracted). Nevertheless, comparison with the component 
field formalism shows that the correct expressions are again reproduced if one takes coin- 
cidence limits after superderivatives. For instance, the correlator {DiX{ti,9i)X{ti,6i)) 
is evaluated by calculating 

{D^x{n,ei)x{T2,92)) = 9^gBi2 - 02^^12, (2.26) 

and then setting T2 = ti. 

It is easily seen that the substitution rule eq.( ^ ) continues to hold, if one defines 
the cycle property solely in terms of the r - indices, irrespectively of what happens to the 
Lorentz indices. For example, the expression 
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^iGsuh £2Qb23£3 hQssih 



would have to be replaced by 



£1^512^2 e2^S23e3 hOsSlh - ei^F12fc2 e2^F23e3 hQpslh 



(2.27) 



(2.28) 



The only novelty is again due to the fact that, in contrast to and Gp, Qb and Qp have 
non- vanishing coincidence limits, 



Qb{t, t) 
Qf{t, t) 



icot{eFT) 
-itan(eFr) 



eFT 



(2.29) 
(2.30) 



As a consequence, we now have also a substitution rule for one-cycles. 



GB{Ti,Ti) GB{Ti,Ti) -GF{Ti,Ti) 



(2.31) 



sin(eFr)cos(eFr) eFT 

This is almost all we need to know for computing one-loop photon scattering ampli- 
tudes, or the corresponding effective action, in a constant overall background field. The 
only further information required at the one-loop level is the change in the path integral 
determinants due to the external field, i.e. the vacuum amplitude in the constant field. 
For spinor QED, this just corresponds to the Euler-Heisenberg Lagrangian, and has, in 



the present formalism, been calculated in |21|. Let us shortly retrace this calculation (the 
fact that the Euler-Heisenberg integrand may be represented as a superdeterminant was 
already noted in f?^). In the scalar QED case, we have to replace 



Vy exp 



dr-y 



Det'p^ 



[A-kT]' 



D_ 

' 2 



by 



Det'p^ 



dl + 2ieFdr = [47rT]-~Det 



1 

P 



2ieFdr~ 



(2.32) 



(2.33) 



(as usual, the prime denotes the absence of the zero mode in a determinant). Application 
of the In det = tr In identity yields 



Det'p' 



I - 2ieFdr^^ 



exp 



exp 



det 2 



--5:^^(-2ie)"tr[F"]TY[ar 



n=l 

oo 



n 



1 Ft 

n — 2 
n even 

"sin(e-FT)' 



eFT 



(2.34) 
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where the are the BernouUi numbers. In the second step eq.( |B.27] ) was used. The 
analogous calculation for the Grassmann path integral yields a factor 



Detp 



det2 



cos(eFT) 



I - 2ieFd~^ 

For spinor QED we therefore find a total overall determinant factor of 

"tan(eFr)" 



[AttT] ~det"2 



eFT 



Expressing tr[F ] in terms of the two invariants of the electromagnetic field, 



(2.35) 



(2.36) 



tr[F 



2ni 



2\" 



- + J (E2 - B2)2 + 4(E • B) 



E^ - B^) + J(E2 - B2)2 + 4(E • B)^ 



(2.37) 



we obtain the standard Schwinger proper-time representation of the (unsubtracted) Euler- 
Heisenberg-Schwinger Lagrangians, 



r(^)\F^ - 1 rdT {eaT){ebT) 
^=^1^^ " (47r)2io sin(ear)sinh(e6r)' 



Ci^)m - _^rdT m^T {eaT){ebT) 

^P^"^^ {AnfJo r3 tan(ear)tanh(e5r) ^ ' 

3. Gauge Boson Loops in External Fields 

In this section we first express the one-loop effective action and also the propagator of 
spin-1 gauge bosons in an arbitrary external Yang-Mills field in terms of a worldline path 
integral. We then will evaluate the action in a covariantly constant background. 

We employ the background gauge fixing technique so that the effective action r[A^] 
becomes a gauge invariant functional of ]73| , 74|. The gauge fixed classical action reads, 
in D dimensions, 

S[a;A] = y d^xF-,iA + a)F-,iA + a) + ^ld^x{D-'[A]a>'^y. (3.1) 

A priori, the background field A'^ is unrelated to the quantum field a^. The kinetic 
operator of the gauge boson fluctuations one obtains as the second functional derivative 
of 5[a, A] with respect to a^, at fixed This leads to the inverse propagator 

P»t = -D^'Dfd^, - 2igF-l (3.2) 
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and the effective action 



r[A] = ^lndet(P) 



In writing down eq. (^]^) we adopted the Feynman gauge a = 1. The covariant derivative 
Dn = dfj_ + igA^j^T"- and the field strength = F^^{T'^Y^ are matrices in the adjoint 
representation of the gauge group, with the generators given by (T"-)^'^ = —if°-^'^. The full 
effective action is obtained by adding the contribution of the Faddeev-Popov ghosts to eq. 
(|3.3| ). The evaluation of the ghost determinant proceeds along the same lines as scalar 
QED and we shall not discuss it here. 

In order to derive a path integral representation of the heat-kernel 

Tr[exp(-rP)] (3.4) 
we first look at a slightly more general problem. We generalize the operator V to 

h^u = -D^5^,u + M^^{x) (3.5) 

where M^^{x) is an arbitrary matrix in color space. In particular, we do not assume 
that the Lorentz trace M^^ is zero. Given M^^, we construct the following one-particle 
Hamilton operator: 

H = {p^ + gA^{x)f- : ^^M,^^, : (3.6) 

The system under consideration has a graded phase-space coordinatized by x^,p^ and 
two sets of anticommuting variables, ip^ and ■0^, which obey canonical anticommutation 
relations: ^ ^ 

For a reason which will become obvious in a moment we have adopted the "anti-Wick" 



ordering in (3.6): all x/j^s are on the right of all ■^'s, e.g. 

: ipf^-ipa ■■ = -^ai>p- (3.8) 

We can represent the commutation relations on a space of wave functions ^{x,ip) de- 
pending on x^ and a set of classical Grassmann variables tp^. The "position" operators 
Xn = Xn, ^11 = tpa act multiplicatively on <I> and the conjugate momenta as derivatives 



= —id^ and V'/t = d/dtp^. Thus the Hamiltonian becomes |77] 



H = _D^ + ^,M,^{x)-^. (3.9) 



The wave functions $ have a decomposition of the form 

D 



p=0 ' 
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This suggests the interpretation of $ as an inhomogeneous differential form on with the 
fermions V'/i playing the role of the differentials dx'^ [76|. The form-degree or, equivalently, 
the fermion number is measured by the operator 

F = i^.i^^ = (3.11) 

We are particularly interested in one- forms: 

$(x,V') =(^^(x)V'^. (3.12) 
The Hamiltonian ( p.S| ) acts on them according to 

{H^){x,^) = {h^,^,)^^. (3.13) 
We see that, when restricted to the one- form sector, the quantum system with the Hamilto- 



nian ( |3.6| ) is equivalent to the one defined by the bosonic matrix Hamiltonian h^^^ 1 76 , 77 ] . 

The euclidean proper time evolution of the wave functions $ is implemented by the 
kernel _ 

K{x2,^2Mxi,^iM) = (a;2,V'2|e-(*2^*^)^|xi,V'i) (3.14) 

which obeys the Schrodinger equation 

d 



^^+Fji^(x,V,T|xo,V'o,0) = (3.15) 

with the initial condition i('(x, V', 0|xo, V'O) 0) = 5{x — xo)6{ip — ipo). It is easy to write 
down a path integral solution to eq. ( p.l5[ ). For the trace of K one obtains 

Tr[e-^^^] = / Vx{t) f Vilj{T)V'4^{T)tvVe' J'o '^'^^ (3.16) 



with 



L = + igA^j_{x)x^ + ipf,[dr6fj,u - M^^]il)y. (3.17) 

We have again periodic boundary conditions for x^[t), and antiperiodic boundary condi- 
tions for iI)^[t). If we use periodic boundary conditions for the fermions we arrive at a 
representation of the Witten index ||7^ rather than the partition function: 

TV[(-1)V^^'^-] = / Vx{t) [ Vi>{T)Vi;{T)tiVe- -fo '^^^ . (3.18) 
Jv Jv 

At this point we have to mention a subtlety which is frequently overlooked but will be 
important later on. If we regard the Hamiltonian (|3.6D as a function of the anticommuting 
c-numbers ^pfJ_ and ipfj_ it is related to the classical Lagrangian (3.17) by a standard Legendre 
transformation. The information about the operator ordering is implicit in the discretiza- 
tion which is used for the definition of the path-integral. Different operator orderings 
correspond to different discretizations. Here we shall adopt the midpoint prescription ||7^ 
for the discretization, because only in this case the familiar path-integral manipulations 
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are allowed Q. It is known ||7|, H] that, at the operatorial level, this is 



equivalent to using the Weyl ordered Hamiltonian H\y. This is the reason why we wrote 
Hw rather than H on the l.h.s. of eqs. ( |3.16| ) and ( 3.18 ). In order to arrive at the relation 
(3.13) we had to assume that the fermion operators in H are "anti-Wick" ordered. Weyl 
ordering amounts to a symmetrization in tp and ip so that 



Hw 



1 



H 



(3.19) 



In the second line of ( 3.19 ) we used (|3.6D and (p?7|). (With respect to and p^, Weyl 
ordering is used throughout.) If we employ (|1|) in (|3l^ ) we obtain the following repre- 
sentation for the partition function of the anti-Wick ordered exponential: 



Ti[e 



-TH} 



Vx{t) / ViIj{t)V^{t) 



X iiV exp 



. dT\L{T) + ]^M^^{x{T))^ 



(3.20) 



Let us now calculate the partition function Tr[exp(— T/i)] which is a generalization of the 
heat-kernel needed in eq. ( |3.3| ). By virtue of eq. ( |3.13D we may write 



Tr[e' 



Tri [e 



-THi 



(3.21) 



where "Tri" denotes the trace in the one-form sector of the theory which contains the 
worldline fermions. In order to perform the projection on the one-form sector we identify 



Mfj_i, with 



M„ 



C6„ 



where C is a real constant. As a consequence, 



H = Ho + CF 



(3.22) 



(3.23) 



with 



Ho = {Pt. + 9Af,{x)f - 2igF,^{x)il,,i!^ (3.24) 

denoting the Hamiltonian which corresponds to the inverse propagator T). The fermion 
number operator F = i^^ip^ is anti-Wick ordered by definition. Its spectrum consists of 
the integers p = 0,1,2, ... D. Note that M^^ = DC, and that because of the antisymmetry 
of F^i, the Hamiltonian Hq has no ordering ambiguity in its fermionic piece. It will prove 



useful to apply eq. ( 3.20 ) not to H directly, but to H — C = Hq + C{F — 1). The result 
reads then 



^[g-CT{F-l)g-THo] 

= exp[-Cr(— - 1)1 / Vx{t) [ Vilj{T)Vi){T)tiVe~ ^0 ^''^ (3.25) 
2 Jp J A 
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with 



L 



(3.26) 



After having performed the path integration in ( |3.25[ ) we shall send C to infinity. While 
this has no effect in the one-form sector, it leads to an exponential suppression factor 
exp[— Cr(p — 1)] in the sectors with fermion numbers p = 2,3, ...D. Hence only the zero 
and the one forms survive the limit C ^ oo. In order to eliminate the contribution from 
the zero forms we insert the projector [1 — (— 1)^]/2 into the trace. It projects on the 
subspace of odd form degrees, and it is easily implemented by combining periodic and 
antiperiodic boundary conditions for ij^^. In this manner we arrive at a representation of 
the partition function of Hq restricted to the one-form sector: 



Tri [e 



lim Tr 

C^oo 

lim exp 

C^oo 



4< 



^CT[F-l)^-THo 



-CT 



D 
~2 



1 



Vx{t) 

)P?/^(r)P^(r)TrPe~ '^''^ 



(3.27) 



Because Tr[exp(— TP)] = Tri[exp(— TZ/q)]) eq. ( 3.27 ) implies for the effective action 



-— lim 

2 C^oo 



X TrP exp 



dT 

T 



exp[-C7r( 



D 



1)] 



dr { ^x^ + igA^x^ + i)^[{dr - C)6^„ - 2igF^y\'iljy 



1 



(3.28) 



Several comments are in order here. The factor exp[— CTL'/2] in ( 3.2g| ) is due to the 
difference between the Weyl and the anti- Wick-ordered Hamiltonian. It is crucial for 
obtaining a finite result in the limit C — > oo. In fact, for D = 4 it converts the prefactor 
exp[+CT] to a decaying exponential exp[— CT] From the point of view of the worldline 
fermions, C plays the role of a mass. Their free Green's function {t2 — Ti)5^y with 

a^(r2 - Ti) =< T2\{dr - C)-Vi > (3.29) 
reads for periodic and antiperiodic boundary conditions, respectively. 



-[G(-r) + e(r)e- 



-CTi 



1 



-CT 



^^(r) = -[e(-r)-G(r)e- 



■CTi 



1 + e- 



-CT 



(3.30) 



We observe that for C — > oo there is an increasingly strong asymmetry between the forward 
and backward propagation in the proper time. Further details of the Green's functions 
(3.30) can be found in appendix B. 

^In ref. |^ the reordering factor was not taken into account and the change of the sign in D = 4 was 
attributed to a difference between Minkowski and Euchdean spacetime which is not correct in our opinion. 
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We mention in passing that there exists another simple method for the projection on 
the one-form sector. We can insert a Kronecker-delta 5^^ into the partition function ( |3.20D 
and exponentiate it in terms of a parameter integral over an angular variable a: 

r2n/T 





T r 







— da e 



iaT 



TV 



^-iaTF^-THo 



(3.31) 



The r.h.s. of ( p. 31 ) can be represented by a path integral which, for the fermions, involves 
antiperiodic boundary conditions only. The corresponding action is similar to the one 
used above but with C replaced by ia. Instead of the limit C ^ oo one has to perform 
the a-integration now. The computational effort is essentially the same in both cases. 

The representation ( |3.28|) of the effective action does not coincide with the one used by 
Strassler |Q. While he uses the same kinetic term in the fermionic worldline Lagrangian, 
he modifies the interaction term according to 



1 



1 



(3.32) 



where 

Xf.{r) = M^) + M^) (3-33) 

As a consequence, the modified Hamiltonian Hq contains terms which change the fermion 
number by 2 units. This means that, during the proper time evolution, the 1-form rep- 
resenting the gauge boson can evolve into a 3-form. However, in the limit C ^ oo the 
substitution ( 3.32|) causes no problems. The reason is that for C ^ oo the energy gap 
between 1- and 3-forms becomes infinite, and therefore a 1-form at r = will remain a 
1-form for alle r > 0. In the modified formalism, Wick contractions of the interaction 
term ( p. 3D involve the 2-point function of X) i-e., ^-^(r) = Q'^'{t) — {—t). From ( |3.30D 
we obtain explicitly 



Qlir) 



sign T 



sign(r) 



sinh.[C(^ - \t\ 
sinh[Cr/2] 

cosh[C7(| - |r| 
cosh[Cr/2] 



(3.34) 



These Green's functions do not coincide with the ones given by Strassler ||9|] 
become effectively equivalent in the limit C 



however, they 
oo. The substitution ( 3.32 ) is motivated 
by the algebraic simplification which it entails in perturbation theory where one inserts 
a sum of plane waves for A'^{x). We shall not do this in the present paper but rather 
calculate the path integral exactly for a covariantly constant field strength. In this case 



the representation ( 3.28 ) is more convenient than the one advocated by Strassler. 

An important building block for higher-loop calculations is the gluon propagator in an 
external Yang-Mills field. It is given by the proper-time integral of the evolution kernel 



( 3.14 ). The latter can be represented by the following path integral with open boundary 
conditions: 



K{x2,il}2,T\xi,ilJi,Qi) 



x{T)=X2 
x{0)=x\ 



Vx{t) 



i/'(0)=i/'i 



PV(t) / 'Dil){T)e 



f'^ drL 
Jo 



(3.35) 
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The Lagrangian L is given by ( 3.17 ) with M^^ = —2igF^y because we do not need the 
C-term in the case of open boundary conditions. The reason is that the Hamiltonian ( |3.9| ) 
preserves the fermion number. Therefore, the kernel K will map a one-form wave function 



of the type ( 3.12| ) onto another one-form. In fact, in the one-form sector, K is represented 
by a Lorentz matrix 



K^y{x2,T\xi,Q) = {x2,fJ.\e 



-TV I 



\Xl,U) 



which is related to the gluon propagator by 

{x2,iJ.\'D^^\xi,u) = / dT K^„{x2,T\xi,0). 
Jo 



(3.36) 



(3.37) 



(The color indices are kept implicit.) It is easy to express the bosonic matrix K^i, in terms 
of the kernel K with fermionic arguments. If one writes 



$(X2,V2,T) = / d^Xid^iJi K{X2,^2,T\XI,^1,0) $(xi,V'i,0) 



(3.38) 



and inserts a wave function of the type ( p. 12 ) at both r = and t = T one finds that 

d 



Kf,i,{x2,T\xi,0) 



dip2,i 



d^tpi K{x2,tp2,T\xi,ipi,0) tpiy. 



(3.39) 



By combining eq. ( 3.37] ) with eqs. (|3.39| ) and ( p. 35 ) we obtain the desired path integral 
representation of the gluon propagator. 

In order to get a better understanding of this representation, let us assume that the 
background is either abelian or quasi-abelian, and that no path-ordering must be 
observed therefore. We may then rewrite ( |3.35 ) according to 



K{x2,iIJ2,T\xi, 11)1,0) 



x{T)=X2 
x(0)=xi 



Vx{t) E:f(V'2,T|v^i,o) 



X exp 



X +igA^x^ 



Here the fermionic integral 

E:F(V'2,r|v^i,o) 



1p(0)=i:i 



X exp 



T 



(3.40) 



(3.41) 



is a functional of the bosonic trajectory x^{t). It can be evaluated exactly [^, 76, ^ and 
has a remarkably simple structure: 



i^F(V'2,r|V'i,o) = ,5(^2 - cS(r)Vi). 



(3.42) 
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In ( 3.42| ) we introduced 



2ig / dTF{x{T)) 



(3.43) 



where the path-ordering is necessary because of the Lorentz matrix structure. Alterna- 
tively, eq. ( p. 421) can be established by noting that Kp solves the Schrodinger equation 
corresponding to ( |3.41| ): 



i^F(V',r|V'i,o) = 0. 



Because = S ^ , it follows that 

and therefore 

Kfj,^{x2,T\xi,0) 



(3.44) 



(3.45) 



x{T)=X2 

X exp 



Vx{t) Pexp 



2ig / dTF{x{T)) 



drl^x^ + igA^x^ 



(3.46) 



Later on we shall use the evolution kernel in the form (3.46). Clearly we could have 
written down this representation without going through the original path integral ( p. 35 ). 
However, in multiloop calculations it will be advantageous if the gluon loops and the 
corresponding propagators are represented in a coherent framework. In fact, if we recall 
that "ip^ amounts to the derivative d/dip^ in the Schrodinger picture, the evolution kernel, 
for any background, may be rewritten in the following very elegant form: 



K^^{x2,T\xi,^) 



x{T)=X2 



Vx{t) / V^{t) / V^{t) 



(3.47) 



In eq. ( 0.47]) , ip{0) and V'(^) are integrated independently. 

Now we evaluate the path integral ( |3.28 ) for the case that the background has a 
covariantly constant field strength. We assume that the gauge field has the form 



Alix) = n''AJx), 



1 



(3.48) 



where is a constant unit vector in color space. The associated field strength = n°'f^u 
with F^^ = d^K - satisfies D^^Fj^^ = 0. Both A^, and F^^, enter eq. (|^ for the 
gauge boson loop as matrices in the adjoint representation. Hence the only nontrivial 
color matrix which enters the path integral is n^T"". If we denote the eigenvalues of this 
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matrix by vi,! = 1,2,..., and evaluate the color trace in the diagonal basis, the integral of 
the path integrand assumes the form 



TrVI{gA^) = J2l{giyiA^) (3.49) 



The path ordering has no effect here. We observe that we are effectively dealing with a 
set of abelian theories whose gauge coupling is given by gi'i. 

While the condition D^F^^ = does not necessarily imply that ^ is constant, we 
further assume that is an x-independent matrix and that the gauge field is of the form 



A^(x) = ^x^F^^ (3.50) 

In the following we keep the diagonalization of the color matrix implicit and continue to 
use the notation TiI{gA^) rather than X^z -^(ff^/A^). We keep in mind, however, that 
and F^y may be treated as pure numbers as far as their color structure is concerned. 
For the field ( 3.4^ ), (|3.5C| ) all path integrals in ( ^.28 ) are Gaussian. We separate the 



constant mode from the integration as before and obtain 
1 /■ ,n dT 



1 



r[A] = -i- I d^xotr / _Deti,[-a2<5^, + 2igF^,dr]-^ lim Y{C) (3.51) 



with 



Y{C) = ^exp[-CT{^-l)]I^DetA[{dr-C)5^u-2igF^,] 

-Betpiidr - C)6^, - 2igF^,]^ (3.52) 

We denote the real eigenvalues of iF^u by f^"\a = 1, and we use the formulas in 

appendix C in order to express Y{C) in terms of these eigenvalues: 

y(c) = iexp[-cr(:^-i)] 



Jn . w.^ r^u 1 IT cosli[^(C + 2g/("))] 
X |Det^[(9. - C) V] n eosh[Cr/2] 

n. l(f^ r^U .1^ sinh[^(g + 2g/("))] l 

-Detp[(a. - c)v] n ,inh[cr/2] — I 

Note that in the case of the fermionic integration with periodic boundary conditions the 
zero mode of dr was not excluded from the determinant and that eq. ( |C.6| ) applies there- 
fore. In ( |3.53| ) we have factored out the free determinants because the method of appendix 
C can yield only ratios of determinants. The normalization factors 

Ca,p = DetA,p[(9, - C) V] = (DetA,p[5, - C]f (3.54) 
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are most easily determined by recalling their operatorial interpretation as the partition 
function and the Witten index of a free Fermi oscillator, respectively: 



BetA[dr-C] =TY[e-^^^^] 
Detp[5, -C]= Tt[{-lfe-^^^^] 



(3.55) 



Here Fw = [tp^ — tpip)/! = F — 1/2 \s the Weyl-ordered fermion number operator (for 
D = 1) with eigenvalues -1/2 and +1/2. As a consequence, 



Ca = (2cosh[Cr/2])^ 
Cp = (2sinh[CT/2])^ 

Taking advantage of Sq,/^"-' = 0, we can rewrite ( 3.53| ) as 

1 



Y{C) = -V 



-1 



+ ma) - 11(1 - r/g. 



(3.56) 



(3.57) 



with r] = exp(— CT) and qa = exp[— 25(T/(°)]. In the limit C oo, i.e., r/ — > 0, the leading 
0(1) terms cancel among the products in ( p. 57 ), and one obtains 



D 

Jim^y(C) = ^exp[-25r/ 

a=l 



(")1 



= tr^cosp^TF] (3.58) 

Here we exploited that if /("^ is an eigenvalue, so is — /'-"^ (tr/, denotes the trace with 
respect to the Lorentz indices.) 

In complete analogy with the scalar case, eq.( ^.35 ), the bosonic determinant in ( |3.51 ) 
gives rise to a factor of 

smigTF)- 



-D/2 



exp 



-— trrln- , 
2 (gTF) 



(3.59) 



Hence our final result for the gauge boson loop becomes 



smigTF)- 



X tr exp 
xtrL cos{2gTF) 



—trrln- , 
2 (gTF) 



(3.60) 



Eq. (3.60) coincides with the result which was found with the help of the traditional 



techniques [^, 84, B5 



By a computation similar to the previous one, but with open boundary conditions for 
the path integral, one can find the gluon propagator in the background (|3.48| ), (|3.50| ). It 
is given by ( 3.37| ) with 



K^i,{x2,T\xi,0) = exp{2igTF) 



x{T)=X2 



x{0)=xi 



'Dx{T)e 



-3[x(-)] 



(3.61) 
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The path integral which remains to be evaluated ist just the one for the scalar propagator. 
Now in the action functional 



(3.62) 



one has to be careful about the surface terms which appear after the integration by parts. 
For open paths they give rise to a non-zero contribution in general. Since S is quadratic in 
Xfj,, the saddle point approximation of the path integral ( |3.61 ) gives the exact answer. We 
can solve it by expanding the integration variable about the classical trajectory connecting 
xi and X2- 

^,(r)=xf-(r)+y^(r) (3.63) 

Here x'^^^^^(r) obeys 



^class 



ligFp^pXy 



•class 



(3.64) 



and it satisfies the boundary conditions x'^^^'^'^(O) = xi and x'^^^^{T) = X2- The fluctuation 
^(t) satisfies correspondingly y(0) = = y{T). Hence the fluctuation determinant is 
almost the same as in the periodic case, the only difference being that there is no zero- 
mode integration in the present case: 



.class" 



K^,y{x2,T\xi,Q) = exp{2igTF)^^ exp(^-S[ 

xDet'p [-d'^5^^ + 2igFf,^dr 

The classical trajectory is easily found: 

class/ X exp(2igFr) - 1 

X (r =xiH rX2-xi). 

ex.p[2igFT) — 1 



-1/2 



(3.65) 



(3.66) 



Its action is entirely due to the surface terms in ( |3.62| ). In Fock-Schwinger gauge centered 
at xi it reads 

^[^class] ^ 1 _ xi)gFcot{gTF){x2 - xi). (3.67) 
Putting everything together we obtain the final result for the propagation kernel: 

Kfj,i/{x2^T\xi^Q) — v**"-^ ; c^PL^'-y-^ j/ii' 

-(x2 - xi)gFcot{gTF){x2 - xi) 



{AttT)- 


-D/2 




■ 1 


X exp 


~4 




■ 1 


X exp 


~2 



tr^ In 



sm{gTF) ' 

iaTF) . 



(3.68) 
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4. A Multiloop Generalization for Quantum Electrodynamics 



We return to the case of quantum electrodynamics, and proceed to the multiloop 
generahzation of the formalism developed in chapter 2. This generalization is constructed 
in strict analogy to the case without an external field [ 36 1 , and we refer the reader to that 
publication for some of the details. 

We first consider scalar electrodynamics at the two-loop level, i.e. a scalar loop with 
an internal photon correction. A photon insertion in the worldloop may, in Feynman 
gauge, be represented in terms of the following current-current interaction term (see e.g. 



1 87, 39, ISH]) inserted into the one- loop path integral. 



e2 r(A) xiTg) • ±in) 



(4.1) 



(A = y — 1). The denominator of this term is again written in the proper-time represen- 
tation, 



r(A) 



A7r^+^{[x{Ta)-x{n)]- 



dT{4:TTT) 2 exp 



(x{Ta) - Xin) 



(4.2) 



It appears then as yet another correction term to the free part of the worldline Lagrangian 
for the scalar loop path integral. It is convenient to rewrite this term in the form 



rT rT 

{x{Ta) - X{n)f = I dTi dT2x{Ti)Bab{Tl,T2)x{T2) 

Jo Jo 



with 



Bab{Tl,T2) = 6{ti - Ta) - S{ti - Tft) d{T2 - To) - 6{t2 - u) 

This allows us to write the new total bosonic kinetic operator as 

92 - 2ieFd - . 



(4.3) 
(4.4) 

(4.5) 



To find the inverse of this operator, we write it as a geometric series. 



(cf - 2ieFd -j;) = (d"^ - 2ieFd) V (^^ - 2ieFd) ^ ^{cl^ - 2ieFd) V • • • , 

(4.6) 

which can be easily summed to yield the following Green's function: 



^(1)/ N r I ^ I ^ [QB{Tl,Ta) 



QB{Tl,n)][QB{Ta,T2) -gB(Tfc,r2)] 
F 9 Cab 



(4.7) 



with the definition 
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Cab = GsiTa^Ta) - GBiTa,n) - GB{n,Ta) + GBin,n) 

cos{eFT) -cos{eFTGBab) 
~ (eFT) sin(eFr) ' 

(4.8) 

Note that this is almost identical with what one would obtain from the ordinary bosonic 

by simply replacing all Gsi/s appearing there by 
the corresponding ^sjj's. The more complicated structure of the denominator is due to 
the fact that the GBij^s are not symmetric anymore, rather we have Qbij = Q'sji^ ^"^^ 
moreover have nonvanishing coincidence limits. The denominator is now in general a 
nontrivial Lorentz matrix, and must be interpreted as a matrix inverse (of course, all 
matrices appearing here commute with each other). 

For the free Gaussian path integral, it is again a simple application of the In det = tr log 
- identity to calculate 



two-loop Green's function |34 



Det'/ 



B 

-d"^ + 2ieFd + — 



Det') 



X Det'p 



Det'; 



TV 



2ieFd~^ 



[47rr]^det 



^ , J2 - 2ieFd 
T 

sin(eFT) 



eFT 



det 



(4.9) 



The generalization from one-loop to two-loop photon amplitude calculations in scalar QED 
thus requires no changes of the formalism itself, but only of the Green's functions used, and 
of the global determinant factor. Of course, in the end three more parameter integrations 
have to be performed. 

As in the case without a background field 
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the whole procedure goes through 
essentially unchanged for the fermion loop, if the superfield formalism is used. As a 
consequence, the supersymmetrization property carries over to the two-loop level, leading 
to a close relationship between the parameter integrals for the same amplitude calculated 
for the scalar and for the fermion loop: They differ only by a replacement of all t/^'s 
by O's, by the additional - integrations, and by the one-loop Grassmann path integral 
factor eq.( 2.35| ). 

The generalization to an arbitrary fixed number of photon insertions is straightfor- 
ward, and leads to formulas for the generalized (super-) Green's functions and (super-) 
determinants identical with the ones given in p^ , ^ up to a replacement of all G^'s (G"s) 
by Qb^s, (^'s). The only point to be mentioned here is that care must now be taken in 
writing the indices of the GBij^s appearing. For instance, the bosonic three-loop Green's 
function must be written 
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1 ^ 



k,l=l 



(4.10) 



The matrix A appearing here is the inverse of the matrix 



Ti - ^(pBaiai — GBaibi — GBb^ai + ^Bfeibi) — ^(^Baia2 " GBaib2 — QBbia2 + ^Bfei 



b2 



— \(GBa2ai — GBa2bi — fefeai + ^Bfebi) ^2 — ^(^£0202 " ^£0262 ~ ^-66202 + ^-86262 

(4.11) 

Ti,T2 denote the proper-time lengths of the two inserted propagators. 

The discussion of the general case of an arbitrary number of scalar (spinor) loops 
interconnected by photon propagators requires no new concepts, and will be deferred to a 



forthcoming review article |9C]. 

While this multiloop construction is done most simply using the Feynman gauge for 
the propagator insertions, other gauges can be implemented as well (the gauge freedom 
has also been discussed in [^]). In an arbitrary covariant gauge, the photon insertion 
term eq.( [4.1[) would read 



^ A-TT 



47r~ -'o 



rf--i 

V 2 



(Xa - XbY 



f-1 



in . fD\Xa-{Xa-Xb){Xa 



Xb) ■ Xb 



{Xa - XbY 



(4.12) 



Here a = 1 corresponds to Feynman gauge, q = to Landau gauge. The integrand may 
also be written as 



r( — 

V 2 



* '^b 



a 



{Xa - Xb) 



fD \ d d \, 



Xb) 



D_ 

' 2 



(4.13) 



This shows that on the worldline gauge transformations correspond to the addition of total 
derivative terms, which is another fact familiar to string theorists (see e.g. [91 1)- This 
form of the photon insertion is also the more practical one for actual calculations. Again 
it carries over to the fermion loop in the superfield formalism mutatis mutandis. 
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5. The Two-Loop Euler-Heisenberg Lagrangian for Scalar QED 



We proceed to the simplest two-loop application of this formalism, which is the two- 
loop generalization of Schwinger's formula for the constant field effective Lagrangian due 
to a scalar loop, eq.( p.38 ). According to the above, we may write the two-loop correction 
to this effective Lagrangian in the form 



(4vr: 



xdet" 



dT 

2^ Jo ~t' 
sin(eFr) 



-■m?Trp- 



oo _ i-T pT 

dT / dTa / dn 

Jo JO 



eFT 



det" 



■'ab 



iVa ■ Vb) 



(5.1) 



We have now a fourfold parameter integral, with T and T representing the scalar and 
photon proper-times, and Ta^b the endpoints of the photon insertion moving around the 
scalar loop. The first determinant factor is identical with the one- loop Euler-Heisenberg- 
Schwinger integrand eq.( |2.35| ), and represents the change of the free path integral deter- 
minant due to the external field; the second one represents its change due to the photon 
insertion. A single Wick contraction needs to be performed on the "left over" numerator of 
the photon insertion, using the modified worldline Green's function eq.( ^]^). This yields 



iVa ■ Vb) 



tr 



Qsab + 



1 {GBaa — G Bab) {Q Bab — QBbb) 



2^ab 



(5.2) 



Care must be taken again with coincidence limits, as the derivatives should not act on 
the variables Tq, Tb explicitly appearing in the two- loop Green's function; again the correct 



rule in calculating {yaVb) is to first differentiate eq.( |4.7|) with respect to Ti,r2, and put 
Ti = Ta,T2 = Tb afterwards. 

After replacing the GbiJ^s and Cab by the explicit expressions given in eqs.( 2.23) and 
eq.( [4.8D , we have already a parameter integral representation for the bare dimensionally 
regularized effective Lagrangian. 

Alternatively one may, in the spirit of the original Bern-Kosower approach, remove Qb 
by a partial integration with respect to Tq or Tb- Using the formula 



and Gsa 



ddet(M) = det(M)tr(dMM"i) 
~^Bba^ one obtains the equivalent parameter integral 



(5.3) 



sin(e-FT) 




eFT 



det 2 



T 



X 2 ^ ti'^Babtr 



Gb 



ab 



T 



2^ab 



+ tr 



{GBaa — GBab){GBab — Gsbb) 



^ 2^ab 



(5.4) 
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For the further evaluation and renormaUzation of this Lagrangian, we will specialize 
the constant field -F to a pure magnetic and to a pure electric field in turns. To facilitate 



comparison with previous calculations |46, 47, 48, 5C], we will moreover switch from 
dimensional regularization to proper-time regularization. This means that henceforth we 
put D = 4, and instead introduce a proper-time UV cutoff Tq later on. 

We begin with a pure magnetic field. The field is taken along the z-axis, so that 
F^'^ = B,F^^ = —B are the only non-vanishing components of the field strength tensor. 
We also introduce the abbreviations z = eBT, and projection matrices 



F = 







10 



0-100 



^ y 



l03 



/ \ 

10 







y 1 y 



Ll2 







10 



10 



y y 



(5.5) 



We may then rewrite the determinant factors eqs.(2.34),(|2.36|) as 



det 2 



det"2 



sin(e-FT) 



eFT 
tan(eFr) 



eFT 



sinh(z) 



tanh(z) 



The Green's functions eq.( 2.19| ),( p.23 ) specialize to 



(5.6) 
(5.7) 



Qb{ti.,T2) 



GB{n,T2) 
Qb{ti,T2) 



T 



cosh(zGsi2) — cosh(2;) 



+ 



2 zsinh(2;) 
T /sinh(zGBi2) 



2z V sinh(2;) 



Gb12 W 



Ll2 



• smh.{zGBi2)^ (cosh.{zGBi2) 1\ .-a 

<J_B12-'-03 H ■ 1 . s -1-12 — . , . s P-r 

sinh(2;j V smh(2;j zj 

T _Lo/^A fcosh(zGBi2)\ zsmh(zGBi2) 

trB12l03 + 2 0i2 ^ . , , . Il2 + 2 „ . . . . iF 



Tsinh(z) 



r ^\ r ^ cosh(zG'Bi2) T ^ 

cosh(2;j 



Tsinh(z) 

sinh(zG'_Bi2 
cosh(z) 



iF 



(5.8) 



In writing Qb have already subtracted its coincidence limit, which is indicated by the 
"bar" . Cah simplifies to 
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Cab = — 2GBafel03 — 2G^Qf,Ii2 , (5.9) 

where we have defined 



We will also use the derivative of this expression, 

sinh(zG'Bafe) 



Similarly we can rewrite 



det 2 



sin(eFT) - 1 
-jp^{T - -Cab] 



1 



tr^Bafetr 



tr Gb 

Gsab 



ah 



T 2 Cab 
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{Caa - Cab)iCab " Cbb) 

\c 



2 '~'ab 



with the abbreviations 



Bab 



sinh(z) 



(5.10) 



(5.11) 



sinh(z) 



77 



3ab 
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ZCOsh{zGBab) 

sinh(z) 

smh{zGBab, 
sinh(z) 



■ sinh(zG'Ba6) -, 

GBabl + \f . 7 

sinh(2;) 



sinh^(zG'Bafe) + cos\i{zGBab) - cosh(z) 



sinh^(z) 



-G%abl 



(5.12) 



7 = (T + Gb. 



ab 



7 



(T + G%ab) 



We rescale to the unit circle, Ta^b = Tua^b-, and use translation invariance in r to set = 0. 
We have then 



GB{Ta,n) = TGBiUa,Ub) =T{Ua- ul) 
GBiTa,n) = GBiUa,Ub) = 1 - 2Ua ■ 



After performance of the T - integration, which is finite and elementary, eq.( |5.4D turns 
into 
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with 



dT 



sinh(z) 



dua A{z, Ua) 



(5.13) 



A 
^3 



L HGBab/G%J A2 A, 

\ ' {GBat - G%^,f ^ iGI,J{GBat " G^^J ^ {GBa,){GBa, " G%J 

4 G%^^,z coi\v{z) - Gsab , 

1 + 2GBahG%ab - ^G%,^Z COth(z) , 
~G^ab ~ '^GBabG%ab ■ 



(5.14) 



Gsab given by eq.(5.10) with T = 1. Here and in the following we often use the 

identity G^^^ = 1 — ^GBab to ehminate GBab in favour of GBab- 

Renormahzation must now be addressed, and will be performed in close analogy to the 
discussion in [^]. The integral in eq.( |5.13| ) suffers from two kinds of divergences: 

1. An overall divergence of the scalar proper-time integral /q°° dT at the lower integra- 
tion limit. 

2. Divergences of the Jq^ dua parameter integral at the points 0, 1 where the endpoints 
of the photon propagator become coincident, Ua = Uh- 

The first one will be removed by one- and two-loop photon wave function renormahzation, 
the second one by one-loop scalar mass renormahzation. As is well known, vertex renor- 
mahzation and scalar self energy renormahzation cancel against each other in this type of 
calculation, and need not be considered. 

By power counting, an overall divergence can exist only for the terms in the effective 
Lagrangian which are of order at most quadratic in the external field B. Expanding the 
integrand of eq.( |5.13| ), K{z,Ua) = ^^^^A{z,Ua), in the variable z, we find 



K{z,Ua) 



12 



GBab G 



Bab 



+ 



1 1 



2 G^ab 



+ 



1 



G 



+ 2 



Bab 



+ 0{z' 



(5.15) 



The complicated singularity appearing here at the point Ua = ui, indicates that this form 
of the parameter integral is not yet optimized for the purpose of renormahzation. In 
particular, it shows a spurious singularity in the coefficient of the induced Maxwell term 
~ z^. This comes not unexpected as the cancellation of sub divergences implied by the 
Ward identity has, in a general gauge, no reason to be manifest at the parameter integral 
level. 

We could improve on this either by switching to Landau gauge, or by performing a 
suitable partial integration on the integrand. The latter procedure is less systematic, but 
easy enough to implement for the simple case at hand: Inspection of the two versions 
we have of this parameter integral, the original one eq.( f>.l\) and the partially integrated 
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one eq.( |5.4| ), shows that we can optimize the integrand by choosing a certain hnear 
combination of both versions, namely 



4cil[^] = ^xeq-(5.1) + ^xeq. 



After integration over T, this leads to another version of eq.( 



„2 ^-oo 
4e / ^-m^T 



with a different integrand 



sinh(z) 7o 



(ilia A'(z,Ua) 



(5.16) 



(5.17) 
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^^{GBahlGBab) , ^^{G Bah I G%ah) 



(<^Ba6)(<^Safe " G|^J {GBab){GBab " G|^j,) 



+ 

3 

A2- 



{GBab — G\ 
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^2 



{G 



Bab 



G 



Bab) 
A' 



tanh(z) 

'-^Bab ~ '-^Bab 



^Bab^Bab "r "-^Bafe 
G\ab + GBahG%ab 



(5.18) 



We have not yet taken into account here the term involving 5ab^ stemming from QBab-, which 
was contained in the integrand of eq.(5.17). This term corresponds, in diagrammatic terms, 
to a tadpole insertion, and could therefore be safely deleted. However, it will be quite 
instructive to keep it and check explicitly that it is taken care of by the renormalization 
procedure. It leads to an integral /g°° ^ which we regulate by introducing an UV cutoff 
for the photon proper-time. 



°° dT _ I 
To f^~%- 

It gives then a further contribution E{Tq) to ^gcait-^]' 

ii;(ro) = -3(4vr)-vi r 

Jo "'0 



(5.19) 



(5.20) 



Expanding the new integrand, K'{z, Ua) 



sinh(2;) 

sinh(2) A-'iz, Ua), in z, we find the simple result 



1 



GBab 



+ + 0(z^ 



(5.21) 



In particular, the absence of a subdivergence for the Maxwell term is now manifest. 
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We delete the irrelevant constant term, and add and subtract the Maxwell term. Defining 



Ko2{z,Ua) = -G-pr ^ 



(5.22) 



the Lagrangian then becomes 



a 



a 



2{47ry Jo 



2{47rf JO 
roo ^j^ 



dT 



dUn 



K'{z,Ua) - Ko2iz,Ua) 



(5.23) 



The second term, which we denote by F, is divergent when integrated over the scalar 
proper-time T. We regulate it by introducing another proper-time cutoff Tq for the scalar 
proper-time integral: 



a 



dT 



(5.24) 



2(4^)^ J2T0 T^' 

(we use 2Tq rather then Tg for easier comparison with |5^ ) . The third term is convergent 
at r = 0, but still has a divergence at Ua = Uh, as it contains negative powers of GBab- 
Expanding the integrand in a Laurent series in GBab-, one finds 



K'(z,Ua) - Ko2iz,Ua) 



Gsab 



+ 



o{g%^,j^ , 



cosh(z) 



sinh(2;) sinh(2;) 



Again the singular part of this expansion is added and subtracted, yielding 



(5.25) 



^SjS] = E{fo)+F{To) 
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^Bab 
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K'{z,Ua) - Ko2{z,Ua) 



GBab 



(5.26) 



The last integral is now completely finite. The third term, which we call G{Tq), is finite 
at T = 0, as f{z) = 0{z'^) by construction. Here we have introduced Tq for the purpose 
of regulating the divergence at Ua = ui,. The Ua ~ integral for this term is then readily 
computed and yields, in the limit Tq — > 0, a contribution 



1 T 



dUa 



1 



GBab 



21n(-^) = -21n(7m2To) + 2\n{-^rr?T). (5.27) 
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We have rewritten this term for reasons which wih become apparent in a moment (ln(7) 
denotes the Euler-Mascheroni constant). Next note that we can relate the function f{z) to 
the scalar one-loop Euler-Heisenberg Lagrangian, eq.( |2.3^ ). If we write this Lagrangian 
for the pure magnetic field case, and subtract the two divergent terms lowest order in z, 
we obtain 



1 



{47ry Jo 



dT 



-m?T 



sinh(z^ 



On the other hand, we can write 



z z^ cosh(z) 
sinh(z) sinh(2;)^ 



d 



1 



z z^ 
sinh(z) 6 



By a partial integration over T, we can therefore reexpress 
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3 rn" dr^_^2r 



(47r)2 7o T2 
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sinh(2;) 6 



(5.30) 



To proceed, we need the value of the one-loop mass displacement in scalar QED, computed 
in the proper-time regularization. This we borrow from [^] 0: 



r 2 3a 2 
dm = — m 

47r 

Using this result, we may rewrite 
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ln(7m^T)/(z). 



(5.32) 



As expected the ^ ~ term introduced by the one-loop mass renormalization cancels the 
tadpole term E{Tq), up to its constant and Maxwell parts. Moreover, the remaining 
divergence of G{Tq) for Tq ^ has been absorbed by 5m'^. 



^Note that this difTers by a sign from 5m^ as used in 
itself, there the corresponding counterterm. 



Here this denotes the mass displacement 
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Putting all pieces together, we can write the complete two-loop approximation to the 
effective Lagrangian in the following way: 
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(5.33) 



We have rewritten this Lagrangian in bare quantities, since up to now we have been work- 
ing in the bare regularized theory. Only mass and photon wave function renormalization 
are required to make this effective Lagrangian finite: 
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e 
B 



ml + 6ml, 
eozl, 



(5.34) 



Here 5ml already been introduced in eq.( |5.31 ), while Z3 is chosen such as to absorb 
the diverging one- and two-loop Maxwell terms in eq. ( 5.33| ). Note that this leaves 
z = cqBoT unaffected. The final answer becomes ^ 
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(5.35) 



^Note added in proof: The constant | multiplying the third term is incorrect, and should be replaced 
by f . This has now been established both by a detailed comparison with and another recalculation 
using dimensional regularization [M. 
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This parameter integral representation is of a similar but simpler structure than the one 
given by Ritus p6[ . 

The corresponding result for the case of a pure electric field is obtained from it by the 
simple substitution (in Minkowski space) 

B -iE. (5.36) 

This makes an important and well-known difference. In the electric field case the T - 
integration acquires new divergences due to the appearance of poles. This leads to an 
imaginary part of the effective action, and to a probability for electron-positron pair 
creation. At the one-loop level, the first such pole becomes significant at the critical field 
strength 

E,r = — (5.37) 

(see eq.( p.38D ). The two- loop contribution to the effective Lagrangian affects also the 
imaginary part and the pair creation probability. A detailed investigation of those correc- 



tions has been undertaken in |47, 45]. 



The calculation for the case of a generic constant field would be only moderately more 
difficult, if one uses the Lorentz frame where the magnetic and electric fields are parallel, 
as in [46|. 



6. The Two-Loop Euler-Heisenberg Lagrangian for Spinor QED 

The corresponding calculation for the spinor loop case proceeds in complete analogy 
when formulated in the superfield formalism. This allows us to immediately write down 



the analogue of eqs.( p3|),( |5.2[ ): 
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with a superfield Wick contraction 
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The notations should be obvious. 

Performing the Grassmann integrations / dOa ^ dOb, and removing GBab by partial in- 
tegration as before, we obtain the equivalent of eq.( 0), 
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(47r 
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(6.3) 



In writing this formula, we have used the symmetry between and to reduce the 
number of terms. The same expressions could have been obtained starting from eq.( 15.4] ) 
and using the generalized one- loop substitution rule. 

Specializing to the pure magnetic case and D = 4, it is then again a matter of simple 
algebra to calculate the traces and T - integrals. After rescaling and setting r;, = 0, the 
result can be written as 



,(2) 
"'spin 
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(47r)" Jto 



dT 
T^' 



-m?T 



tanh(2;) 



dUaB{z,Ua), 



(6.4) 



with 



1 [Gsab - G^^f^f G%^f^{GBab - G^^i,) GBabiGsab " G^^J j ' 

Bi = ^1 - 4ztanh(z)G^„fe , 
B2 = ^2 + 8ztanh(z)G|j„j,-3 , 
^3 = ^3 -4ztanh(2;)G|j„fe + 3 . 

(6.5) 

Comparison with an earlier field theory calculation performed by Dittrich and one of the 



authors |5C] shows that the integrand of eq.( p^ ) allows for a direct identification with 
its counterpart there, as given in eqs. (7. 21), (7. 22). This requires nothing more than a 
rotation to Minkowskian proper-time, T — > is, a transformation of variables from Ua to 
V := Gsabi and the use of trigonometric identities. In particular, our quantities Gsab-, G%^^ 
then identify with the quantities a, h there. 

The renormalization of this Lagrangian has, for the spinor-loop case, been carried 
through in detail in that work. We will therefore not repeat this analysis here, and just 
give the final result for the renormalized two-loop contribution to the Euler-Heisenberg 
Lagrangian ^ : 

Note added in proof: The constant —10 m ultip lying the third term is incorrect, and should be replaced 



4 



by —18 MM (compare the footnote before eq.(5.35)) 
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(6.7) 



For a study of the strong field Hmit of this Lagrangian see again |5C]. 

The first exact calculation of this two-loop Lagrangian for the fermion loop case is again 



due to Ritus |49|, who also used proper-time methods to arrive at a certain two-parameter 
integral. 

The integral representation given above is equivalent to the one given by Ritus, but 



simpler. In |5C] it was obtained by convoluting a free photon propagator with the polariza- 
tion tensor of a fermion in a constant magnetic field. The essential part of this calculation 
consists of deriving a compact integral representation for the polarization tensor. To this 
end, complicated expressions involving Dirac traces and momentum integrals have to be 
evaluated. In the string-inspired case, no analogous manipulations are needed, and the 
computational effort for doing the parameter integrals which it introduces instead is much 
smaller. 

Moreover, when applied to spinor QED, the method of the present paper yields the 
corresponding result for scalar QED with almost no further effort. This would not be the 
case for the standard field theory techniques. 

Concerning the physical relevance of this type of calculation, let us mention the exper- 
iment PVLAS in preparation at Legnaro, Italy, which is an optical experiment designed 
to yield the first experimental measurement of the Euler-Heisenberg Lagrangian p^, 94]. 



It is conceivable that the technology used there may even allow for the measurement of 



the two- loop correction in the near future |94, 95 1. 
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7. The 2-Loop QED f3 Functions Revisited 



Finally, let us remark that the method we have employed in this paper for the calcu- 
lation of the full two-loop Euler-Heisenberg Lagrangians also improves on the calculation 
of the two- loop QED /3 - functions as it had been presented in [^] . For the extraction of 
the P - function coefficients one needs only to calculate the induced Maxwell terms. Up 
to the contributions from one-loop mass renormalization, the correct two-loop scalar 196 



and spinor |97| QED coefficients can thus be read off from the expansions (see eq.( p.21\} 
eq.(0)) 

K'{z,Ua) = --^+3z2 + 0(z4)^ 
<~rBab 
1 2 

L{z,ua) = --^+2^2 + 0(^4). 

(7.1) 

Comparing with [36| we see that the use of the generalized Green's functions Qb,Gf has 



saved us two integrations: The same formulas eq.( |B.2| ) which there had been employed for 
executing the integrations over the points of interaction ti,T2 with the external field, have 
now entered already at the level of the construction of those Green's functions. Of course, 
for the f3 - function calculation all terms of order higher than 0{F'^) are irrelevant, so that 
one could then as well use the truncations of those Green's functions given in eq.( p. 241 ). 
Moreover, one would choose an external field with the property F'^ ~ I. 

Note that in the fermion loop case a sub divergence- free integrand was obtained pro- 
ceeding directly from the partially integrated version eq.( [6.31) . This fact, which had 
already been noticed in [^, is not accidental, and can be understood by an analysis of 
the quadratic divergences. In the scalar QED case, there are three possible sources of 
quadratic divergences for the induced Maxwell term: 

1. The contact term containing 5ab- 

2. The leading order term ~ — in the expansion of the main term (see e.g. 



eq-O). 



Bab 



3. The explicit I/Tq appearing in the one-loop mass displacement eq.( p. 31 ). 



The last one should cancel the other two in the renormalization procedure, if those are 
regulated by the same UV cutoff Tq for the photon proper-time, and this was verified in 
various versions of this calculation. In the spinor QED case the fermion propagator has 
no quadratic divergence (this is, of course, manifest in the first order formalism, while 
in the second order formalism there are various diagrams contributing to the one-loop 
fermion self energy, and the absence of a quadratic divergence is due to a cancellation 
among them). The third term is thus missing, and the other two have to cancel among 
themselves. In particular, the completely partially integrated version of the integrand has 
no 5ab ~ term any more, and consequently the second term must also be absent. However 
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the ^ - expansion of the main contribution to the Maxweh term is, if one does this 
calculation in = 4, always of the form shown in eq.( ^.15 ), 



A B ^ 
+ ^ + C 



-'-'Bab '-'Bab J 

with coefficients A, B, C. In the partially integrated version first consideration of the 
quadratic subdivergence allows one to conclude that A = 0, and then consideration of the 
logarithmic subdivergence that i? = 0. 

Note that this argument does not apply to the scalar QED case, nor does it to spinor 
QED in dimensional renormalization, due to the suppression of quadratic divergences by 
that scheme. In both cases one would have only one constraint equation for the two 
coefficients A and B appearing in the partially integrated integrand, and indeed they 
turn out to be nonzero in both cases. In the present formalism, the fermion QED two- 
loop (3 ~ function calculation thus becomes simpler when performed not in dimensional 
regularization, but in some four-dimensional scheme such as proper-time or Pauli-Villars 
regularization. 

The reader may rightfully ask why we have gone to such lengths in analyzing this 
2-loop calculation, which is easy to do by modern standards even in field theory. We 
find this cancellation mechanism interesting in view of some facts known about the three- 
loop fermion QED /3 - function 1 95 , |100|] . Apart from the well-known cancellation 



of transcendentals occuring between diagrams in the calculation of the quenched (one 



fermion loop) contribution to this (3 - function [pS] , IOC], which takes place in any scheme 
and gauge, even more spectacular cancellations were found in [^] where this calculation 
was performed in four dimensions, Pauli-Villars regularization, and Feynman gauge. In 
that calculation all contributions from nonplanar diagrams happened to cancel out exactly. 
A recalculation of this coefficient in the present formalism is currently being undertaken 

m. 



8. Discussion 

In this paper, we have extended previous work of two of the authors on the multiloop 
generalization of the string-inspired technique to the case of quantum electrodynamics in a 
constant external field. The resulting formalism has been tested on a recalculation of the 
two-loop corrections to the Euler-Heisenberg Lagrangian for quantum electrodynamics. 
Several advantages of this calculus over standard field theory methods have been pointed 
out. In particular, it treats the scalar and spinor loop cases on the same footing, so that 
the scalar loop results are always obtained as a byproduct of the corresponding spinor loop 
calculations. More technically, our parameter integrals are written in a form convenient 
for partial integrations. In particular, the integrands are functions well-defined on the 
circle, so that boundary terms do not appear. The usefulness of this property has been 
demonstrated in the renormalization of the scalar QED two-loop Lagrangian. 

An application to a recalculation of the one-loop QED photon splitting amplitude has 



been given separately |102]. 



We have derived a path integral representation of the gluon loop, and used it for a 
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recalculation of what is the closest analogon to the one-loop Euler-Heisenberg Lagragian 
in Yang-Mills theory. More significantly, the analysis of chapter 3 should be viewed as 
a first step towards an extension of the worldline technique to multiloop calculations in 
nonablian gauge theory. Our derivation of this path integral is entirely non-heuristic, 
and thus guaranteed to reproduce the correct one-loop off-shell amplitudes for Yang-Mills 
theory. From our experience with quantum electrodynamics, this property alone makes us 
optimistic about the existence of a multiloop generalization of the worldline method for 
the Yang-Mills case. We hope to have more to say about this in the future. 
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Appendix A: Path Integral Representation of the Electron Propagator 



Different from the case of the closed loop fermionic worldline Lagrangian eq. ( , for 
the Dirac propagator besides the einbein field gauged to T one has to introduce a gravitino 
field X ill the world line action. This can be gauged to a constant, but not to zero. For 



massive fermions a further field ^5{t) coupling to mass has to be introduced |13, |14|, 19] 
(its super symmetric partner X5 is not needed for the gauge coupling, but essential for the 
worldline implementation of Yukawa couplings ||2^ ) . 

Litegration over x in the path integral leads to a factor (— ^V'AtX^+mV's) corresponding 
to the numerator of the Dirac propagator. This can be demonstrated nicely for the free 



propagator |55, 56, 27 1 in the coherent state formalism. This "holomorphic representation" 
represents operators in terms of their (fermionic) Wick symbols. 

In the case with background interaction considered here we prefer a different approach 
in which fermionic operators are represented by their Weyl symbols. This formalism is 
manifestly covariant and, contrary to the holomorphic representation, it treats propagators 
in external fields and one-loop effective actions on the same footing. From a canonical 
point of view we are dealing with the following algebra of hermitian operators tp^^ and -05: 



In terms of euclidean Dirac matrices, it can be represented by 

i " 1 

V'm = ■^^S^A'' ^5 = (A.2) 

The Weyl symbol map "symb" establishes a linear one-to-one map between operators 
and functions of the anticommuting c-numbers and ^5. In particular, symh{ipf^) = 
and symb(V'5) = Cs- The inverse symbol map associates the Weyl-ordered (totally 
antisymmetrized) operator product to strings of ^'s: 

symb"^(^^.^,,...^p) = {ipf,i>u--i'p}weji (A.3) 

For example, 

symb"^(^pCj.) = -{i^i^i^v - ^ui^p) = i^pi^u - -V (A-4) 

and SYmh{il)f^ipiy) = ^^^y + (See ref. [103,104] for further details.) 

Let us consider the Dirac propagator in the background of an arbitrary abelian gauge 
field and let us write down a path integral representation for its kernel (bosonic vari- 
ables)/symbol (fermionic variables) 

GD''^^=(x2,xi;0 = symb[(x2](/) + m)-Vi)](6 (A.5) 

with p = 'JpDp = 2iip^'il)^D^. After having integrated out the auxiliary fields x s-iid X5 it 
reads [13,27,104], up to an overall constant: 
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G^''-''%X2,xi;0 oc r dTe-"'''^ (A.6) 
Jo 

x{T)=X2 r r 

Vx / Vi)^ / PV's 

x(0)=xi JV(0)+i/'(T)=25 -'V'5{0)+V.5(T)=2C5 



1 1 

/ dT{--'4)^{T)Xf,{T) + rnV'sC'r)} exp[-S'B - 5"^ - ^5 
i JO ^ 



The action consists of the following pieces: 

/o 



-Sb = dT[-xl + ieAf^{x)Xf, 



1 . 1 

Sf = / drl-^p^Tp^ - ieF^^{x)Tpf,il)^] + 
Jo ^ 2 

dr^V5^5 + ^V'5(T)V5(0) (A.7) 

Note the surface terms in Sp and S5. They are needed in order to correctly reproduce the 
equations of motion [104]. The factor y Jq drj...} in eq. ( |A.6| ) stems from the integration 
over the world line gravitino field x- It is important to realize that the terms inside the 
curly brackets are actually independent of r, and that we may replace 



The r-independence of the expectation value of iIj^x^ is a consequence of the supersym- 
metry of the action. In fact, Sb + Sp is invariant under 

5x^ = -27]t/j^ 

5ip^ = r]Xfj, (A. 9) 

with a constant Grassmann parameter tj. If, instead, a time-dependent parameter is used 
in ( A. 9 ), the action changes according to 

S{Sb + Sf) = r dT r]{T) -^(^^i^) (A.IO) 

JO ClT 

Obviously V'/^^;^ is the conserved Noether charge related to the supersymmetry ( |A.9D . If we 
apply a localized supersymmetry transformation to the path integral J DxDip exp(— S b — 
Sf) and observe that the measure is invariant we obtain the Ward identity 

^{i;^{r)x^{T))=0 (A.ll) 



Eq. ( A.lip together with a similar argument for ^5 justifies the replacement ( [A.Sl ). 

Using ( |A.8| ) in (|A.6|) , the insertion —■^ip/iX^ + rm/j^ is evaluated at the final point of 
the trajectory, t = T. Hence it may be pulled in front of the path integral, then acting as 
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a (differential/matrix) operator on the wave function which was time-evolved by the path 
integral. If we are dealing with a phase-space path integral of the type 



x{T)=X2 
x{0)=x\ 



Vx / Vp exp{ / dT{ipx — H)} 



we know that 



J VxVp xiT) exp{...} = X2 j VxVp exp{...} 

/ VxVp p{T) exp{...} = —i— — / VxVp exp{...} 
J OXo J 



(A.12) 

(A.13) 
(A.14) 



By rewriting eq. ( A.6 ) in hamiltonian form, it is easy to see that in the case at hand x^{T) 
corresponds to the operator —2iD^{x2) acting from the left. With an analogous reasoning 
for the fermions this leads us to the following representation for the Dirac propagator: 



Dirac 



G 



Here we defined 
and 



{x2,xi) oc [■4^f,iD^{x2) +m'ip5] 

dTe-"'^ symb-^fi^D''^^^(x2,r|xi,0;C/.) hi^5,T) 



(A.15) 



x{T)=X2 
x{0)=xi 



Vx 



Vip^ e~ 



^(0)+V(T)=2^ 



(A.16) 



(A.17) 



Eq. ([Al5|) was obtained from ( |A.6| ) by applying the inverse symbol map. As for the 
fermionic degrees of freedom, G^^'^'^'^{xi, X2) is an operator now, i.e., a matrix acting on 
spinor indices. 

Up to this point, no assumption about the gauge field Afj^{x) has been made. From now 
on we consider fields with F^j^ = const. In this case the path integral (|A.16|) factorizes: 



K^''^^X2,T\xi,0;^, 



The bosonic piece 



KBix2,T\xi,0) 



Kb{x2,T\xi,0) If{^^,,T) 

x{T)=X2 



Vx e 



-Sb 



(A.18) 



(A.19) 



' x{0)=xi 

is the same as in the spin-0 or spin-1 case. Its evaluation is described in detail in section 
3. The result is given by (3.68) with g replaced by e, and with the factor exp[2igTF]^jy 
omitted. What remains to be done is to calculate the fermionic contribution 



PV'm exp[-i^^(r)V';.(0)] 
V'(o)+v(r)=2e 2 



X exp 



(A.20) 
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Since the ■0-integral is Gaussian, the saddle point method will yield its exact value. We 
decompose the integration variable according to 



where is a solution of the classical equation of motion, 



(A.21) 



(A.22) 



subject to the boundary condition ip'r}'^^ (0) + 'ip'r^'^^{T) = 2^^. The fluctuation field 99^ 



class / 



satisfies antiperiodic boundary conditions. Using (|A.21| ) and ( |A.22 ) in ( A.2C| ) we obtain 

If{C,,T) = exp[-i<-^(r)V'^i-^(0)] 

= exp[-iv'f'^(r)V^i-^(0)] DetA[5,5^,-2ieF^,]V2 (A.23) 

In analogy with section 3, the determinant in (A.23) is given by 2^detL[cos(eFr)]. The 
solution to (|A.22 ) which satisfies the correct boundary conditions reads 



exp(2ieFr) \ 
l + exp(2ieFr)j^/'' 



(A.24) 



Inserting this into ( A.23 ) leads us to the final result for Ip: 

If{^^,T) = 2^/2 detL[cos(eFr)]V2 
X exp[i^^tan(eFr)^,,^i,] 



(A.25) 



Using the same method we can show that I5 equals an unimportant constant which we 
shall drop. Thus, because 



symh-^K^''^'') = Kb symb-\lF), 



(A.26) 



our last task is to find out which is the operator corresponding to the symbol ( A. 251) . We 
shall see that 

symb[exp(zerF^^'0^^^^)](^) 

= detL[cos(eFr)]i/2 exp^^ tan(eFr)^^e..] (A.27) 

In order to prove ( A.27| ), we transform F^^, to block-diagonal form We assume D even. 



F, 



diag 



B^' 



B2\ 
-B2O ' 



(A.28) 
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Since the V'^^'s pertaining to different blocks are mutually anticommuting, we may prove 
( A.27| ) for each block separately. Focusing on the first one, it is convenient to define 



Si2 = i(-0iV'2 
As this operator is Weyl-ordered, we have 



symb 



■'12 



(A.29) 
(A.30) 



Because (^12)^ = 1, it follows that 
exp [zeTF^,, -0^7/^1,] 



exp[eSiTSi2] 

cosh(ei?iT) — S12 sinh(ei3ir) 



(A.31) 



and therefore 



symb[exp(ieTF^,,-(/'^^/ii,)] (^) 

= cosh(eSir)[l - 2i^i6 tanh(e5ir)] 
= det L [cos(eFr)] exp [i^^ tan(eFr) 



(A.32) 



In the last line of ( |A.32D we used that the eigenvalues of the first block are ibii?i and that 



^\ = = ^2- Repeating the same argument for the other blocks establishes eq. ( A.27| ) 



Upon inserting ( |A.26| ) with (|A.25| ) and (|A.27] ) into (|A.15|) , we obtain the well-known 
result for the euclidean Dirac propagator in a constant background field [104,105]: 



G°*'-^'(xi,X2) = [-7^Z)^(x2)+m] / dT{A^T) 



-D/2g-m2T 



• exp 



• exp 



• exp 



'^(3;2 — 2;i)eF cot(eTF)(a;2 — 2:1) 

1 , sinfeTF) 
--trL In— i 

2 (eTF 



(A.33) 



In ( A.33| ) we used the representation ( A.2| ) for the "^'s, and we dropped an overall factor 
of 75 which is produced by the path integral, but is not included in the standard definition 
of G^^^^. The expression for the bosonic contribution Kb was taken from section 3. It 
applies to the Fock-Schwinger gauge centered at xi. In the general case, Kb contains an 
extra phase factor exp(— ze /^^^ dx^A^). We also note that the scalar propagator (— Z)^)^^ 
is obtained from ( A.33| ) by simply deleting the operator [—p + m] and the last exponential 
involving the 7-matrices. 

It is remarkable that the above calculation of the propagator is almost identical to 
the calculation of the one-loop effective action, the only difference being the boundary 
condition of the fermionic path integral. For the propagator we need (T) -|- -i/^^j (0) = 2^^, 
whereby the variables give rise to its 7-matrix structure. The effective action, on the 
other hand, is a scalar quantity, and it is obtained from the same path integral with = 0. 
Giving a non-zero value to amounts to creating a fermion line by "opening" a loop. 
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Appendix B: Derivation of Worldline Green's Functions 

The worldline Green's functions appearing in this paper are kernels of certain integral 
operators, acting in the real Hilbert space of periodic or antiperiodic functions defined on 
an interval of length T. We denote by Hp the full space of periodic functions, by Hp 
the same space with the constant mode exempted, and by Ha the space of antiperiodic 
functions. The ordinary derivative acting on those functions is correpondingly denoted by 
dp, dp or Ba- With those definitions, we can write our Green's functions as 



Qb{ti,T2) = 2(ti I (dp^ - 2iFdpy^ I rs), 

GFin,T2) = 2(ti I [Oa - 2iF)"' I T2), 
G^{ti,T2) = in I {dp-Cy^ I T2), 

Q^{ri,T2) = in I {dA-Cy^ I rs). 

(B.l) 

(in this appendix we absorb the coupling constant e into the external field -F). Note that 
Q'^ is, up to a conventional factor of 2, formally identical with Qp under the replacement 
C 2iF. 

Qb and Qp are easy to construct using the following representation of the integral 
kernels for inverse derivatives on the unit circle [36| 



(u I ap-" I u') = -—,BJ\ u-u' Dsign'^fn - u') 
m 

{u\dA-''\u) = —^En-l{\u-u'\)s\gn^{u-u') 

2[n — 1)\ 



(B.2) 



Here Bn{En) denotes the n-th Bernoulli (Euler) polynomial. Those formulas are valid for 
\ u — u' \ < 1. For instance, the computation of Qp proceeds as follows: 



Qb{ui,U2) = 2(ni I (ap2 - liFdp) ^ I U2) 

00 

= 2^ (2iF)"('Ui|ap'("+2)|^2) 

n=0 

°0 /o-TTiN" — 2 



(2iFr-^sign"(ni-^2) p „ 

2 > 1 Bn(\ Ui - U2 \) 

n=2 

1 signfiii — n2)e^*-^("i~"2) signfui — U2) „ 

H 7-p; -Did ""1 - ""2 



g2jFsign(«i-U2) _ I IP iVI i z 1; 2^2 

1 



2F2 V sinF 



^-iFGB,2 + iFGpu - 1 



(B.3) 



43 



In the next-to-last step we used the generating identity for the Bernoulh polynomials, 



= J2Bn{x)- (B.4) 
1 n] 



n=0 

This is Gb as given in eq.( p.l9 ) up to a simple rescaling. The computation of Qp proceeds 
in a completely analogous way. 

This method does not work for the determination of Q^,, as negative powers of Bp are 
not even well-defined in the presence of the zero mode. 

In the following, we will calculate Qj, _^ in a different, more "physical" way, which 
corresponds to the usual construction of the Feynman propagator in field theory. 

In order to determine G^{t), say, we employ the following set of basis functions over 
the circle with circumference T: 



/„(T)=T-i/2exp[i^(n+i)r], nGZ (B.5) 



They satisfy 







T 

dTf*{T)frri{T) = 5„ 



E /n(r2)/:(ri) = E 5{T2-n-mT) (B.6) 



m=— CO 



and /„(t -|- r) = —fn{T)- In this basis, the Green's function ( |3.29[ ) becomes 
By introducing an auxiliary integration in the form {t = ti — T2) 

/OO 1 °° / 9-77- 1 \ 

-00 T V T 2 Jiu-C 

and using Poisson's resummation formula, the Green's function assumes the suggestive 

00 

GA{r)= E (-ira£(r + nr) (B.9) 



form 1 106] 



with 



We verify that 



1^-^ (B.IO) 

^00 2,11 lijj — C 



[dr - C]Gg{T) = 5{t) (B.ll) 

00 

[dr-C]gg{T) = E (-im^ + nT) (B.12) 



44 



which shows that is a Green's function on the infinitely extended real line, while 
is defined on the circle. The integral ( [B.10| ) yields for C > 

a£(r) = -e(-r)e^- 



Hence, from ([B.9| 



gC^r) = -e^^ E {-ire{-T - nT)e 



nCT 



(B.13) 



(B.14) 



For r G (0,T) only the terms n = — oo,...,— 1 contribute to the sum in ( |B.14D , while 
for r G (— T, 0) a nonzero contribution is obtained for n = — oo,...,0. Summing up the 
geometric series in either case and combining the results we obtain the expression given 
in eq. ( 0.30 ). It is valid for —T < r < +T. Using a basis of periodic functions the same 
arguments lead to Qp as stated in (p.SOl). Note that in the limit of a large period T 



lim g^p{T) 



(B.15) 



as it should be. For C — > 0, both and G'^ have a well-defined limit: 



9\{r) - 



-G(-t) 



sign [T) 



(B.16) 



The periodic Green's function Qp blows up in this limit because dp^ does not exist in 
presence of the constant mode. It is important to keep in mind that Gp is defined in such 
a way that it includes the zero mode of dr- 

In the perturbative evaluation of the spin-1 path integral one has to deal with traces 
over chains of propagators of the form 



Because 



(B.17) 



(B.18) 



(n- 1)! \dC, 

it is sufficient to know p{C). The subtle point which we would like to mention here is 
that strictly speaking the sum defining a\, say, 



4(c) = E 



i{2'n /T){n + \) - C 



(B.19) 



does not converge as it stands, and it is meaningless without a prescription of how to 
regularize it. The usual strategy is to combine terms for positive and negative values of 
n, and to replace ( B.19 ) by the convergent series 



4(C) 



-2CE 

n=0 

— tanh , 
2 V 2 



I) 



/CT 



(B.20) 
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It is important to realize that this definition imphes a weh-defined prescription for the 



treatment of the functions in p at r = 0. In fact, 



a\{c)= r dr g1{T-T) = Tg^{0), 
Jo 



and by combining eqs. ( |B.2C| ) and ( B.2l| ) we deduce that we must set 

1 



hm e(r) = hm Q{-t) 

r\0 t\0 Z 



With ( |B.20|) we obtain 



1 



T\" / d 
dx 



{n-iy. \ 2 

The analogous relation in the periodic case is 



tanh(x^ 



x=CT/2 



1 /T 



(n - 1)! V2 



dx 



n-l 



— — coth(x' 



x=CT/2 



(B.21) 



(B.22) 



(B.23) 



(B.24) 



if the zero mode of dr is included in the trace (B.17), and 



(n-l)! V2 



dx 



n-l 



{coth(2;) — X ^} 



x=CT/2 



(B.25) 



if the zero mode is omitted. For C sufficiently small one finds the power series expansions 



1 ^ (2 ' - l)B2k j,2k(j2k~n 



B2k rp2k(j2k-n 



(B.26) 



a\ and a p have well defined limits for C ^ 0: 



(2" - 1 En-l n I \ 

-1 = —1 [n even) 



n 



2 {n-l) 



T" (n even) 



(B.27) 



Those limits vanish for n odd. This brings us, of course, back to eqs.( [B.2| 
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Appendix C: Worldline Determinants 



In this appendix we collect a few results about the determinants which arise in the compu- 
tation of the spin-1 wordline path integral. To start with, we consider the operator d-r + uj 
where a; is a real constant, and dr acts on periodic and antiperiodic functions of period 
T, respectively. Its spectrum reads i{27r/T)n in the former and i{2iT/T){n + 1/2) in the 
latter case, n £ Z. Ratios of determinants of the form 

DetA[dr + io] ^ i(f )(n + |)+^ 

are defined by the prescription that terms with positive and negative values of n should 
be combined so as to obtain the manifestly convergent product 



BetA[dr+uj] _ ^ 



1 + 



LvTy 1 



27r J (n + 1/2)2 



cosh f ^) (C.2) 



In the periodic case we omit the zero mode from the definition of the determinants and 
find likewise 

Det'pjdr+u;] ^ sinh(a;r/2) 

Det'p[dr] {ujT/2) ^ ' 

(compare eqs.( |2.35| ), ( |2.35| )). Next we look at the matrix differential operator 

{dr - C)5^y + Vl^^, /X, = 1, D. 

Here is a constant matrix. We assume that it can be diagonalized and has eigenvalues 
uj. For antiperiodic boundary conditions we obtain from (|C.2D 

T>(itA[{dr-C)5^,^ + n^^] cosh[|(C - a; 



DeiA[{dr - C) V] y cosh[Cr/2] ' ^^'^^ 

The product extends over the spectrum of Vt. The corresponding formula for periodic 
boundary conditions, with the zero mode removed, reads 

Det'p[{dr-C)5^,] lic-a;'-^ sinh[CT/2] ' ^'^ 

For C 7^ we can reinstate the zero mode of dr- In this case ( p. 5]) is replaced by 



Detp[(9^ - C) V + n^u] _ sinh[|(C - uj)] 

Bet p[{dr - C)S^u] y sinh[Cr/2] ' ^'^ 

In this paper we use the above determinants for an exact evaluation of the worldline 
path integral in the background of a constant field Ffj^i^. For more complicated field con- 
figurations only a perturbative calculation of the path integral is possible in general. It 
is based upon the Green's functions Gap which were discussed in appendix B. It can be 
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checked that the determinants given above are consistent with the perturbative expansion. 
In perturbation theory, the l.h.s. of eq. ( p.4| ), for instance, is interpreted as a power series 
in u): 



[] DetA[l + ^Qa\ = n TV hi[l + wC?^] 
= Jl exp 



(C.7) 



n=l 



In the last hne of ( |C.7| ) we have used ( B.17| ). By virtue of eq. ( |B.23 ) one can sum up the 
perturbation series in closed form: 



n=l 



n 



^ 2"n! \dxj ^ 

n=l ^ ' 

In cosh(3;) 



exp 
In 



LuT d 
~Ydx 
cosh[x — ujT/2] 



x=CT/2 
x=C'T/2 



cosh(x) 



x=CT/2 



(C.8) 



With (p.8D inserted into (C.7) we reproduce precisely the r.h.s. of eq. 
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